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Plasma Wakefield Effects On High-Current 
Relativistic Electron Beam Transport In The Ion-Focused Regime 

J.D. Miller, R.F. Schneider, D.J. Weidman<'‘\ and H.S. Uhm 
Naval Surface Warfare Center, Silver Spring, MD 20903-5000 

and 

K.T. Nguyen 

Mission Research Corporation, Newington, VA 22122 

Abstract 

Modulation of the beam current has been observed during ion-focused regime (IFR) 
transport of a high-power relativistic electron beam propagating through a low-density 
b'’Ci;ground plasma. Injecting a high-current, high-energy electron beam into an IFR chanre’ 
immersed in a background plasma induces plasma oscillations. These background plasma 
oscillations, induced by the risetime portion of the beam ejecting plasma electrons from the 
vicinity of the beam into the background plasma, give rise to a modulated axial electric field. 
This field travels with the beam leading to beam energy and current oscillations. In the 
experiment, a 1.7-MeV, 1-kA, risetime-sharpened electron beam is propagated on a KrF excimer 
laser-produced IFR channel in TMA gas, which is immersed in a low-density plasma-filled 
transport tube. We present experimental measurements and computer simulations demonstrating 
modulation of this high-current relativistic electron beam near the low-density background 
plasma frequency. 

PACS numbers: 41.80.Ee, 52.40.Mj 
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I. INTRODUCTION 


Plasma wakefields are well known in the high-energy physics community. The plasma 
Wakefield accelerator (PWFA)‘ utilizes the longitudinal electric fields in plasma waves excited 
by a short duration bunched relativistic electron beam as a source of high accelerating gradients. 
In the PWFA, intense electron bunches are utilized as a driver beam to excite wakefields in a 
high-density plasma. These fields are then used to accelerate a low-current trailing beam.^’^ We 
report here experimental measurements on high-current (kA level) relativistic electron beams 
whose transport is strongly influenced by the self-induced wakefields. 


A high-current relativistic electron beam injected into a preionized plasma channel causes 
channel electrons to be expelled by the electrostatic force generated by the head of the beam and 
permits stable electron beam transport.*** If the plasma channel is immersed in a low-density 
background plasma, as illustrated in Fig. 1, both channel electrons and plasma electrons out to 
the charge neutralization radius are ejected by the beam head, provided the combined channel 
and background plasma density is less than the beam density. The charge neutralization radius 
is the point at which the total enclosed channel and background plasma ion charge is equal to 
the beam charge. The charge neutralization radius can be easily shown to be 
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where r^ is the beam radius, r,. is the channel radius, fc=ni/nb is the space charge neutralization 
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fraction and is the ratio of the channel ion density to the electron beam density, and g=np/nb is 
the ratio of the plasma density outside the channel to the electron beam density. The background 
electrons beyond the charge neutralization radius are not expelled, but they are perturbed by the 
beam head and begin to oscillate at a frequency near that of the electron plasma frequency of 
the low-density background plasma, a)p=(npeVeome)‘^. These radial plasma oscillations produce 
an electrostatic wakefield with electric field components in the radial and axial directions. The 
longitudinal component of this field, E,sE<^in(a)pt), traveling with the beam, causes accelerating 
and decelerating forces on beam electrons leading to beam energy and current modulation. These 
effects were first observed for high-current beams in computer simulations® and were seen to 
lead to eventual beam disruption. Recently, direct experimental measurements have been made 
which confirm the existence of this phenomena.*” 

The intensity of the energy modulation and the resulting degree of current modulation 
through the beam pulse is determined by the strength of the wakefield. For an electron beam 
propagating through a uniform background plasma contained within a grounded conductor of 
radius r„, the amplitude of the axial electric field has been determined, via computer simulation 
using the two-and-one-half dimensional PIC code MAGIC**, to scale as 

= 0.18 n^*/^ 4 [F/cm], 

where np is in units of cm'^, and Ib is the beam current in kA. This scaling is valid for beam 
currents in the range Ib < 2 kA and beam rise times less than the background plasma period. 
This scaling indicates that the higher the plasma density, the stronger the beam modulation for 
the same propagation distance. Numerical simulations also indicate that the amplitude of the 
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plasma wakefield is reduced from that shown in Eq. (2) as the beam rise-time is increased. For 
long-range propagation, the initial beam rise-time is eroded away after some distance and the 
wakefield amplitude is again increased. For the present experimental configuration, the 
simulations have indicated that it will be necessary to sharpen the beam rise-time prior to 
injecting the beam into the wakefield chamber. 


II. EXERIMENTAL ARRANGEMENT 


In the experiment, a 1.7-MeV, 6-kA, 30-ns electron beam produced by a Febetron 705 
is passed through a beam risetime sharpener and matched onto a laser produced IFR channel in 
trimethylamine (TMA) gas. Figure 2 is a schematic diagram of the experiment. The electron 
beam is generated from a 2.5-cm-diam velvet cathode with a 2.2 cm anode-cathode gap. The 
electron beam is extracted into a beam risetime sharpening cell through a 6 ^tm aluminized mylar 
anode foil. The risetime sharpener cell uses a single magnetic lens to preferentially focus the 
high energy portion of the beam at an aperture further downstream. This downstream aperture 
is a 3-cm-diam graphite aperture plate covered by a second aluminized mylar foil. The beam 
exiting the risetime sharpening cell is matched onto a KrF excimer laser ionized IFR channel in 
TMA. The TMA gas flowed through the transport chamber continuously to provide a stable gas 
pressure in the range 0.2 - 0.8 mTorr. These TMA pressures result in space-charge neutraliza¬ 
tion fractions sufficient to efficiently propagate the electron beam in the absence of a background 
plasma. The TMA gas pressure is low enough that beam-induced ionization from the REB is 
negligible. The beam risetime sharpener cell and the diode region are maintained at a lower 
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pressure, on the order of lO"* Torr, to prevent diode shorting or gas focusing effects. The laser 
is fired to ionize the TMA typically 500 ns before the voltage was applied to the diode. The 
resulting laser produced channel is centered in a 0.5-ni-diam, 3.6-m-long plasma filled transport 
chamber. The low-density background plasma filling the transport chamber is generated by a hot 
filament discharge in the low pressure TMA gas. 

The beam current is monitored prior to injection into the plasma filled transport chamber 
and after exiting the chamber by resistive wall current monitors.*^ Uncalibrated, single-turn B- 
dot loops*’ located axially along the transport chamber monitor the evolution of the modulation 
of the current. Diametrically opposing B-dot loops could also be configured to detect transverse 
beam centroid motion.*’ It should be noted that both current monitors are sensitive to the net 
current flowing in the transport chamber. A scintillator-photodiode observed the x-ray signal 
emitted when the beam electrons were deflected to the wall at the end of the transport chamber 
to avoid striking the laser input window. 

A. Beam Risetime Sharpening 

The electron beam generated by the Febetron 705 nominally has a pulse shape that is 
roughly triangular both in voltage and current. The diode voltage is monitored with a capacitive 
voltage divider. The diode current is monitored with a Rogowski coil located immediately 
downstream of the anode foil. Typical beam voltage and diode current waveforms are shown in 
Fig. 3. The electron beam current risetime of 15-ns is several times the plasma oscillation period 
of the low-density background plasma filling the transport chamber. Ideally, a zero risetime 
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beam would result in maximum wakefleld excitation.*'* A 2 ^ro risetime beam is difficult, if not 
impossible, to achieve in practice however. The risetime can be sharpened sufficiently to 
enhance the wakefield excitation to an observable level. This is desirable in the laboratory 
experiment due to the finite propagation distance available. 

A simple method for shaping the output current pulse of a relativistic electron beam in 
vacuum has been implemented to sharpen the risetime of the Febetron generated beam.*^ The 
technique utilizes a magnetic lens with a magnitude of about 1-kG to focus the beam. Passing 
beam electrons through the magnetic lens causes them to focus at different axial locations 
downstream from the lens depending upon their energy. The focal point of the beam current 
peak (corresponding to maximum energy) is then located furthest downstream. By placing a 
small aperture near the beam peak focal point, we let the portion of the beam near the peak pass 
through, while intercepting most of the early portion of the pulse which has spread out before 
it reaches the aperture. Similarly, the later portion of the pulse fall time will also be intercepted. 
This creates an electron beam with sharper rise and fall times resulting in a more uniform 
current profile. This technique results in a 1-kA electron beam exiting the risetime sharpening 
cell with a 5-ns risetime and a 13-ns flattop. 

B. Laser Ionized Channels in TMA 

The 3-cm-diam laser ionized channel was formed by two-step photoionization of TMA 
using a 750-mJ, 248-nm, 30-ns KrF excimer laser. Measurements were performed*® to 
characterize the laser photoionization of: (a) initially neutral TMA gas, and (b) TMA gas 


6 



1 


partially ionized by a heated filament DC discharge operated in both a pulsed (ms timescale) and 
steady-state mode. 

Well-defined ionization volumes have been demonstrated to be quite efficient using 
multiphoton ionization of organic molecules. In this experiment we use a high-power KrF 
excimer laser to ionize a channel in trimethylamine [(CH 3 ) 3 N] gas (TMA). Previous 
measurements*’ have indicated that TMA has a high two-photon ionization cross-section at the 
KrF laser wavelength of 248 nm. In our application it is necessary to determine accurate 
parameters to allow the degree of ionization to be calculated for arbitrary laser properties. A 
simple model that agrees well with previous experimental data for fractional ionization yields 
on the order of or less, is that the fractional ionization (the ratio of plasma density Up to the 
neutral gas density n,,) depends quadratically on the laser fluence.^** This relationship takes the 
form 


Hi = CF^ ( 3 ) 

where F is the laser fluence (mJ/cm^) and C is the absolute ionization yield coefficient. 
Generally C is defined as C=of/hj'T, where o; is the two-photon ionization coefficient, hp is the 
photon energy, and t is the laser pulse half-width. Previous measurements have determined the 
value for C to be in the range 1.52x10^ to 1.83x10-® cmVmJ^ for laser pulse widths between 
13.2 and 5.65 ns respectively for KrF laser ionization of TMA.^® 

The situation of particular interest to our experiment is whether this relationship is 
appropriate for laser photoionization of partially ionized TMA. The photoionization process may 
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be effected due to the presence of a low-density DC discharge. In addition, the presence of such 
a discharge may lead to fragmentation of the large organic molecules and a corresponding 
reduction in the photoionization efficiency. 

The experimental setup used to measure the laser parameters was similar to that shown 
in Fig. 2. The beam size at the output window of the laser is approximately 1 cm x 2.7 cm. The 
measured beam divergence was found to be very close to the specified value of 0.2 mrad. The 
laser pulse width is monitored with a DC biased pin diode and has a FWHM of approximately 
25 ns. The laser beam is passed through a prism to expand the cross-section to about 3.1 cm 
square, and then through two adjustable apertures in series to limit the beam cross-section in the 
test cell. The laser intensity was then assumed to be spatially uniform in the test region. The 
laser energy per pulse is measured using a Scientech calorimeter. The TMA pressure was 
measured with ionization gauges calibrated to a Baratron capacitance manometer. The test cell 
was initially evacuated to about 4x10^ Ton* and TMA was allowed to flow continuously through 
the chamber. 

The degree of ionization produced in the chamber by the laser pulse is determined by 
collecting the charge obtained from a known laser volume on biased electrodes inserted into the 
test cell.^* The results of the ionization yield measurements performed in neutral TMA are shown 
in Fig. 4, which plots the fractional ionization as a function of the laser fluence. The 
measurements show the expected F^ dependence over the range of the data. The absolute 
ionization yield coefficient is extracted by least squares fitting the experimental data to Eq. (3). 
This fit is shown as the solid line in Fig. 4. From this fit we have determined C to be 4.27x1(1’ 
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cmVmJ^. This result is in reasonable agreement with the previous measurements^® and supports 
the use of Eq. (3) to determine the plasma density for fractional ionization yields up to 10 ^ 


The presence of the DC steady-state discharge influences this result considerably. The 
measurements made in the presence of the low-density discharge plasma utilized microwave 
resonator probe techniques^* and simple DC biased Langmuir probes^^ collecting electron 
saturation current. The partially ionized plasma background was produced with a heated filament 
discharge operating in both a DC steady-state mode and a DC pulsed (ms timescale) mode and 
will be described in the next section. Figure 5 demonstrates the reduced photoionization yield 
of the TMA measured by the microwave resonator probe in the DC steady-state discharge. The 
fractional ionization is seen to be reduced by about a factor of two in the low-current discharges. 
The photoionization yield then seems to be less sensitive to further increases in discharge 
current. The solid lines are curves fitted to an exponential decay of the form n=noe‘'% 
characteristic of plasma loss due to diffusion. This effect is also observed with the DC biased 
Langmuir probe. 

This result has been found to be due primarily to the fragmentation of the TMA 
molecules in the low-density steady-state plasma discharge. Measurements of the fragmentation 
of the TMA in the discharges was made with the quadrupole gas analyzer. A typical mass 
spectra of the constituents present in the neutral TMA gas is shown in Fig. 6(a). The TMA line 
(59 amu) is most prominent with the existence of a few impurities at much lower partial 
pressures. The presence of the DC steady-state discharge modifies the mass spectrum 
considerably as shown in Fig. 6(b). Most apparent is the reduction of the TMA line and the 
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increase in the impurity lines to about the same level. The increase in impurities is accounted 
for primarily by H, H2, N, CH3, and multiples of these. These are readily identifiable as 
fragments of TMA. Addition of argon (40 amu) as a buffer gas for the discharge did not aid in 

reducing the fragmentation as shown in Figs. 6(c) and (d). 

The presence of the DC steady-state discharge resulted in an unacceptibly large reduction 
in the photoionization yield of the TMA. For the relativistic electron beam transport experiments 
the plasma channel parameters must be independent of the low-density background plasma 
conditions. To reduce the fragmentation of the TMA, the low-density background plasma 
discharge was operated in a pulsed DC mode. In this case the high-voltage discharge was pulsed 
for times less than 1 ms. This resulted in significantly reducing the degree of TMA fragmenta¬ 
tion as measured with the DC biased Langmuir probes*®. 

C. Large Volume Plasma Formation 

To fill the relatively large transport chamber uniformly with low-density (10* - 10^ cm *) 
background plasma, we have utilized a technique based upon direct ionization of the background 
gas by low-energy electrons emitted from heated tungsten filaments.^* Heated filament discharges 
are a very basic and well-established technique for plasma generation. An annular arrangement 
of the filaments produces a relatively uniform radial plasma density profile. Application of a 
weak axial confining magnetic field (on the order of 5 G) ensures axial plasma uniformity and 
provides plasma confinement. This field is sufficienfly small so as not to influence the electron 
beam dynamics. The discharge is pulsed for durations on the order of 1-ms to reduce the level 
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of TMA fragmentation before the laser is fired.** The Langmuir probe data indicates that the 
plasma is formed within the first 200 and behaves as a steady-state for the remainder of the 
voltage pulse. Properties of the plasma are diagnosed using Langmuir probe techniques and a 
microwave resonator probe^* located at diagnostic ports as shown in Fig. 2. 

Peak plasma densities as a function of the discharge current are shown in Fig. 7. The 
relationship between the plasma density and the discharge current is found to be approximately 
linear for these conditions. Note that the error bars on the Langmuir probe data in this and 
subsequent figures are indicative of relative errors only. The microwave resonator probe 
provides the more accurate absolute plasma density measurements. Electron temperature 
measurements made with the Langmuir probe ranged from 1 to 3 eV. 

A combination of Langmuir probe geometries was used to provide axial, radial, and 
azimuthal probe adjustability. Data for the axial, radial, and azimuthal plasma density profiles 
were taken individually over several minutes while the discharge was operating in the DC 
steady-state mode using argon as the fill gas. The gross features of these plasma density profiles 
were later confirmed with the discharge operating in TMA. The radial and axial plasma density 
profiles under a variety of discharge conditions are shown in Figs 8(a) and 8(b) respectively. 
The radial plasma density profile is observed to be quite uniform inside the annular filament 
location before falling off close to the walls of the chamber. Similarly, the axial plasma density 
profile taken along half the length of the chamber shows good uniformity over 1.4 m before 
increasing somewhat near the end of the chamber where the filaments are located. Since the 
geometry of the chamber is symmetric on both ends, the midplane has been defined as the zero 
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point on this plot and measurements taken along one end only. The plasma density is adjustable 
in the range 10* to 5x10’ cm *. 


III. EXPERIMENTAL RESULTS 


Initial experimental measurements have demonstrated very well defined current 
oscillations impressed upon the electron beam after traversing the 3.6-m-long transport chamber. 
These oscillations are in the range 150 to 300 MHz and are very close to the electron plasma 
frequency of the low density background plasma. In the absence of the background plasma, the 
most efficient electron beam transport ocurred for a space-charge neutralization fraction of 
feS»0.9 as determined by the laser channel ionization measurements. This is indicative of a fairly 
high transverse beam temperature and is common to electron beams generated and propagated 
using these techniques.*^ 

The evolution of the current oscillations is characterized in the typical waveforms shown 
in Fig. 9. Figure 9(a) shows the response of the resistive wall current monitors in the absence 
of the background plasma. The currents are shown prior to the entrance of the beam into the 
transport chamber (upper waveform) and after exiting the chamber (lower waveform). 
Examination of the transported current waveform indicates a steepening of the current risetime 
with an associated spike at the front end of the pulse. This is due to a combination of inductive 
beam erosion and the fact that the beam voltage pulse is nearly sinusoidal and varies throughout 
the current pulse.** Figures 2(b)-(d) illustrate the development of current oscillations with 
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increasing background plasma density. The plasma densities have been inferred from DC biased 
(+40 V) Langmuir probe data without any correction for finite sheath effects. The space charge 
neutralization fraction at the peak beam current for all the data presented here was fc~0.9, 
giving a neutralization radius, r„, ranging from 3.0 cm (Fig. 9(b)) to 2.4 cm (Fig. 9(d)). For 
these experimental parameters, n/nb assumes values in the range 0 < np/nb ^ 0.06, indicating 
that the low density background plasma does not contribute significantly to the space charge 
neutralization fraction of the channel. The data in Figs. 9(c) and (d) indicate that the current 
becomes severely modulated after only 3.6 m of transport in the plasma background. 
Measurements performed with the resistive wall current monitors and the B-dot loops configured 
to detect transverse beam motion observed no discemable indication of beam motion transverse 
to the direction of propagation. Increasing the plasma density higher than that shown in Fig. 9(d) 
resulted in loss of transport efficiency and eventual beam disruption as indicated by the current 
monitors and the x-ray signal picked up by the scintillator-photodiode. Although a useful 
indicator for gross beam dynamics, the response time of the scintillator was not fast enough to 
respond to the beam current modulation. 

The frequency of oscillation is more readily determined from the unintegrated B-dot loop 
data. Figure 10 shows the unintegrated B-dot loop response from the probe located 0.9 m 
downstream from the entrance to the plasma-filled transport chamber. The current oscillations 
are quite clear and indicate an increasing frequency with increasing plasma density. The 
dependence of the background plasma density on the oscillation frequency obtained from the 
unintegrated B-dot loop data is shown in Fig. 11. The solid line represents a least squares fit of 
the data to a rip'” dependence. This functional dependence indicates that the experimentally 
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observed current oscillation frequency scales quite well with the low density background electron 
plasma frequency. Comparison of the experimentally observed oscillation frequencies with the 
electron plasma frequency of the low-density background plasma inferred from the DC biased 
Langmuir probe data, indicates that the experimentally observed oscillation frequency is between 
75 to 80% that of the electron plasma frequency calculated from the Langmuir probe data. A 
quantitative comparison such as this however is difficult to state with a high degree of 
confidence due to the uncertainty in the plasma density inferred from the Langmuir Probe data 
(which may be as high as a factor of 2). 


IV. COMPARISON WITH THEORY 


A. Axisymmetric Plasma Wakefield Mode 


In this section we theoretically investigate the current profile of the relativistic electron 
beam upon exiting the plasma filled transport chamber. We assume that the axial electric field 
of the Wakefield has a sinusoidal waveform with axial wavenumber k. Due to inductive erosion 
at the beam head, the beam front velocity j8,c is slightly less than the initial beam velocity ^t,c. 
Labeling to for the beam segment which enters the chamber at time t=to, we obtain the axial 
momentum change Ap(z,t)^* 


eE^ 


t4>{zj) = ^ p [^cos^O-^APcfo) - cos(0-A:Apcr)j 


(4) 


when the beam segment to propagates to z at time t. Here the normalized velocity difference 
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Aj3=i8b-|8f, Eq is the maximum electric field of the wakefield given by Eq. (2), and the angle 
parameter 6 is defined by 0=k/?bCto. In our experiment, the combination of beam current and 
relatively short propagation distance make the inductive erosion rate rather small. In the absence 
of significant erosion, i.e., the momentum change in Eq. (4) is simplified to 


AP(z,^o) = " 




sinO. 


( 5 ) 


In reality, the velocity of the beam segment to is continuously changing during the 
propagation and this velocity is expressed as ^c=p/ytn, where the relativistic mass ratio y is 
defined as 7 =(l+p^/mV) ‘'^, and the axial momentum p and the mass ratio y are functions of 
the position z and time t. The velocity change during the propagation is then related to the 
momentum change Ap by 6/3 = Ap/ 7 ^mc. The normalized axial velocity ^ of the beam segment 
to is therefore expressed as /3=/3b+6i3, when it arrives at z. By approximating the mass ratio 7 
by 7bj where 7 b is the initial beam energy, we can obtain the velocity from Eq. (5), 




1 - 


eE^ 




sinO 


( 6 ) 


of the segment to when it arrives at z. The approximation of 7 by 7 b is valid only for 7 b> > 1. 
Note that the velocity of the segment to is j 8 c=/ 8 bC at z=0 and that ^b(to) is ^ function of to in 
general. Thus, the average velocity of the segment to over the range from z=0 to z=z is given 
by 


(p)c = P j,c (1 - K sin0) 


( 7 ) 


where the parameter k is defined by 
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Note that the initial beam current Ib(to) enters continuously through the chamber entrance 
located at z=0 and time t=to. When this beam segment arrives at z=z in time t, this beam 
segment is stretched by the factor dt/dto- Therefore, the beam current of the segment to at z is 
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proportional to d0/d<^. Making use of the normalization factor 


1 } ^ = 1 ( 12 ) 

71 J 1 + fl cosx 

we obtain the beam current I(tQ) 

M. = ^ (13) 

IJto) l+i«cose’ 

assuming that there is no beam loss during the propagation. In reality however, we can expect 
beam current loss and approximate this with a fractional propagation factor a, which is in the 
range 0<a< 1. In obtaining Eq. (13), we have assumed that the parameter k is less than unity, 
which is valid for a limited propagation distance. 


We have made several assumptions in the preceding analysis. First, beam-head erosion 
is negligibly small. For the beam parameters in the experiment, this is a reasonable approxima¬ 
tion due to transient effects resulting from the relatively short propagation length. Second, 
although beam segments can be stacked up to each other, individual beam segments do not cross. 
Third, the beam segment gains or loses energy according to a predetermined sinusoidal 
waveform. In reality, the energy-gain mechanism of a beam segment may be much more 
complicated than a simple sinusoidal wave pattern. This simple theoretical description does 
however provide a basic understanding of the beam current profile at the end of the chamber. 


For comparison, we consider the case for the background plasma density np—10 cm , 
which corresponds to the plasma frequency Wp=l.78x10* rad/s. The wakefield frequency is 
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expressed o=^(jip, where the crossover parameter q is typically 0.65. Thus, the 

Wakefield frequency is 1.43x10^ rad/s. We also assume that the beam energy is about 1.1- 
MeV, which represents an average value of the beam energy throughout the beam pulse. The 
relativistic mass ratio is 7b=3.2. The beam current at the injection point is 500 A and the space- 
charge neutralization fraction is fe=0.5. The beam risetime is assumed to be t,=5 ns. These 
conditions are a good approximation to the experimental realization described previously. The 
maximum electric field Eq of the wakefield wave is not sensitively dependent upon the beam 
radius for beam risetimes longer than one plasma oscillation period.*'* Making use of these 
conditions, we find Eo=2.8 kV/cm from Eq. (2). Estimating the chamber length to be about 3.5 
m, we find that /c=3xl0'^. Note that the value of the parameter k is very small. The parameter 
X is given by x=0.53. 

Shown in Fig. 12 are plots of the input and output current ratio versus 6 

obtained from Eq. (13) for x=0.53, and a=l and a=0.5. Here a=l represents the case of no 
beam electron loss and ex =0.5 means that one-half of the beam electrons are lost during 
propagation. As expected, the output beam current has very sharp peaks at d=(2n+ 1)t, where 
n=0,l,2,... is an integer. Note that this simple theoretical estimation of the output current 
profile is very similar to the experimentally observed data. 

It is interesting to note that for long distance propagation, in the case where the 
parameter x is larger than unity, the solution indicates that the current I(to) has two cusped peaks 
for each period of the wakefield oscillation. Symptoms of this behavior have been observed in 
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the particle simulation studies^® of long distance beam propagation and also appear to be evident 
in the present experiment and simulations. 

B. Numerical Simulations 

Numerical simulations have been performed to predict the behavior of the electron beam 
in the plasma filled chamber using the fully electromagnetic 2-1/2 dimensional particle-in-cell 
code MAGIC.“ This code was developed for simulating axisymmetric plasma phenomena such 
as the case at hand. The beam and channel properties in the simulation correspond quite closely 
to the experimental conditions described previously. The simulation parameters are listed in 
Table 1. The beam initial energy had an isosceles triangle profile with a base of 640 kV and a 
peak of 1.7 MV. The transport chamber was filled with a uniform diffuse plasma with a density 
that could be varied, although unless otherwise indicated, the simulation results shown here will 
be at a background plasma density of 2x10^ cm'^. 

Shown in Fig. 13 is an electron trajectory plot consisting of both beam electrons and 
background plasma electrons. The electron beam is propagating from left to right in each frame. 
In frame (a), the trajectory of the electrons are shown 4.8 ns after injection into the plasma 
chamber. As expected, the beam head is expanding radially since there is no electrostatic 
focusing force. Further back from the beam head, as the channel and background plasma 
electrons are ejected from the vicinity of the beam the beam is able to propagate in a focused 
manner. Figure 13(b) shows the evolution of this behavior at 9.6 ns. In this frame, oscillations 
of the background plasma about the charge neutralization radius are clearly apparent. By 28.8 
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ns into the pulse (Fig. 13(d)), the beam is actually being chopped up into beamlets. 

The ability of the diffuse plasma Wakefields to strongly modulate the beam energy is 
clearly shown in the beam axial momentum plot in Figure 14(a). The wakefield is the result of 
the diffuse plasma oscillations seen in Fig. 13(a)-(d). The axial length of the wakefield shown 
in Fig. 14 corresponds very favorably with the plasma wavelength of 75 cm expected for a 
diffuse background plasma density of 2x10’ cm l This can also be observed in the axial electric 
field plot at the corresponding time (t=19 ns) as shown in Fig. 14(b). One particular feature 
which can be noted from Fig. 14(b) is the growing amplitude of the wakefield as a function of 
axial distance. We believe this is due to the fact that as the beam energy is modulated by the 
wakefield, the beam bunches as it propagates downstream. The modulated beam current 
resonantly pumps the diffuse plasma, enhancing the wakefield amplitude, which in turn 
accelerates the beam bunching process. 

The axial evolution of this bunching process on the electron beam current in the transport 
chamber is shown in Fig. 15(a)-(d). By the time the electron beam reaches the end of the 3.6 
m transport chamber, the beam current is quite strongly bunched. This result compares very well 
with the experimental data shown in Fig. 9(d). The number of peaks (bunches) observed in the 
simulation of Fig. 15(d) is in excellent agreement with the experimental observations. It is 
interesting to note that the bunching frequency observed in the simulation is about 320 MHz, 
or 20% below the expected plasma frequency of 400 MHz for np=2.0xl0’ cm'^. Also notable 
in the frames of Fig. 15 is the generation of modulated plasma currents after the end of the beam 
current pulse. This effect is also quite apparent in the experimental data of Fig. 9. It is 
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worthwhile to point out that the relatively larger peaks of the current in the simulations are due 
to the fact that the current is integrated over the whole cross-section of the wakefield chamber. 
In the experiments, the current was measured with a resistive wall current monitor located at the 
edge of a smaller (15 cm diameter) tremspoit tube a short distance downstream of the wakefield 
chamber. 

A more direct comparison between the simulation results and the experimental data is 
shown in Fig. 16. The beam current profiles in Fig. 16(a)-(d) correspond closely to the 
experimental conditions shown in Fig. 9(a)-(d). The simulation clearly shows the ability of the 
diffuse plasma to strongly modulate the beam current. The simulation result compares very well 
to the experimental data. Figure 16(a) even shows the current spike at the front of the beam due 
to the combination of inductive erosion and the triangular beam voltage profile. 


V. SUMMARY 

We have observed very well-defined current oscillations impressed upon a high-current 
relativistic electron beam propagating in an IFR channel immersed in a low-density background 
plasma. Severe effects on beam current transport have been observed for g parameter values in 
the range 0.03 to 0.06 after only 3.6 m of beam transport. Initial results indicate the frequency 
of current oscillation to be near the plasma frequency of the low-density background plasma, 
consistent with the predictions of plasma wakefield theory. Numerical simulations using the 2- 
1/2 dimensional particle-in-cell code MAGIC are in very good agreement with the experimental 
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measurements. 
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TABLE CAPTION 


Table 1. MAGIC simulation parameters. 




FIGURE CAPTIONS 


Figure 1. Pictorial representation of an electron beam propagating on an IFR channel 
immersed in a low-density background plasma. 

Figure 2. Schematic of the plasma wakefield experiment. 

Figure 3. Diode voltage and current waveforms. 

Figure 4. Measurement of photoionization yield for TMA: o parallel plate ionization 
detector, □ microwave resonator probe. 

Figure 5. Photoionization yield for TMA in a steady-state, low-density dc discharge plasma. 

Laser fluence is 38.1 mJ/cm^ for the hollow symbols, and 36.9 mJ/cm^ for the 
solid symbols: o no background plasma, v np=2.3xl0® cm'*, ° np=5.6xl0* 
cm'*, • no background plasma, ▼ np=2.3xl0* cm *. 

Figure 6. Mass spectra of a steady-state dc discharge in TMA; (a) no plasma discharge, 
TMA pressure 6.8x10'* Torr, (b) 1.0 A plasma discharge current, 6.8x10'* Torr, 
(c) no plasma discharge, TMA pressure 6.8x10'* Torr with 1.1x10'^ Torr Ar, (d) 
1.0 A plasma discharge, TMA pressure 6.8x10* Torr with 1.1x10^ Torr Ar. 

Figure 7. Plasma density as a function of discharge current. Discharge voltage = -70 V, 
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confining magnetic field = 5 G. Hollow symbols are argon at SxlO"* Torr, filled 
symbols are TMA at 6.6x10"* Torr. Circles represent Langmuir probe data, 
squares represent microwave resonator probe data. 


Figure 8. Normalized (a) radial plasma density profiles in argon; discharge voltage is -70 
V; o 2.0 A discharge current, 5 G magnetic field, 2x10^ Torr, z=0 m, 
normalized to 4.4x10* cm'*, ° 1.25 A discharge current, 7.5 G magnetic field, 
2x10-^ Torr, z=0 m, normalized to 5.0x10* cm'*, a 2.0 A discharge current, 5 
G magnetic field, 2x10^ Torr, z=l.l m, normalized to 5.3x10* cm'*, and (b) on 
axis axial plasma density profiles in argon: discharge voltage = -70 V, magnetic 
field = 5 G, o 2.0 A discharge current, 2x10"* Torr, normalized to 4.0x10* cm *, 
° 3.0 A discharge current, 5x10"* Torr, normalized to 6.0x10* cm'*. 

Figure 9. Current waveforms from the resistive wall current monitors; (a) np=0, (b) 
np=1.0x10’ cm'*, (c) np= 1.5x10’ cm'*, (d) np=2.0xl0’ cm*. Top (bottom) 
waveform: current entering (exiting) the transport chamber. 

Figure 10. Unintegrated B-dot loop response 0.9 m downstream from the entrance to the 
plasma filled transport chamber; (a) np=0, (b) np=6.3xl0* cm'*, (c) np=1.0x10’ 
cm'*, (d) np=1.5x10’ cm *. 

Figure 11. Oscillation frequency from the unintegrated B-dot loop data as a function of the 
background plasma density. The solid line is a least squares fit to Up*'*. 
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Figure 12. 


Figure 13. 


Figure 14. 


Figure 15. 


Figure 16. 


Plots of the output to input beam current ratio versus 6 obtained from Eq. (13). 

Plasma and beam electron trajectories from the MAGIC simulation for parameters 
as listed in Table 1. 

Electron beam axial momentum (a) and the axial electric field (b) at t= 19 ns from 
injection into the plasma chamber. 

Beam current waveforms from the MAGIC simulation at (a) z=0 m, (b) z=1.2 
m, (c) z=2.4 m, and (d) z=3.6 m downstream from the entrance to the plasma 
filled transport chamber for np=2.0xl0’ cm'^. 

Beam current,waveforms from the MAGIC simulation for (a) no background 
plasma, (b) 11^=1.0x10’ cm■^ (c) np=1.5xl0’ cm•^ and (d) np=2.0xl0’ cm■^ 3.6 
m downstream from the entrance to the plasma filled transport chamber. 
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SYMBOL 


PARAMETER 


VALUE 


Ij, beam current 1.0 kA 


r^o initial beam radius 1.5 cm 


V beam voltage 1.7 MV 


beam risetime 5 ns 


relativistic mass factor 4.5 


kTj^ initial transverse temperature 1.657 keV 


channel neutralization factor 0.8 


r^. channel radius 1.5 cm 


background plasma density 0-2x10’ cm’^ 


r^ Wakefield chamber radius 25 cm 


L^^, Wakefield chamber length 3.6 m 


































































fractional ionization (n^/n 











OOOVO^fMO OOOVO^<NO 

^OOOOO T-OOOOO 


(qje) 3inss3jd |e{)jed (qje) eanssajd lej^jed 















plasma density (cm-3) ; 








normalized plasma density normalized plasma density 










current (500 A/div) current (500 A/div) 



10 ns/div 


10 ns/div 






unintegrated B-dot signal unintegrated B-dot signal 
(0.5 V/div) (0.5 V/div) 



10 ns/div 


10 ns/div 










2.0 












































-0.4 



-1.50 I-'-1-^-1 -3.0 





































-0.25 



time 
































